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We  present  a  Raman  scattering  study  of  wurtzite  ZnO  over  a  temperature  range  from  80  to  750  K.  Second- 
order  Raman  features  are  interpreted  in  the  light  of  recent  ab  initio  phonon  density  of  states  calculations.  The 
temperature  dependence  of  the  Raman  intensities  allows  the  assignment  of  difference  modes  to  be  made 
unambiguously.  Some  weak,  shatp  Raman  peaks  are  detected  whose  temperature  dependence  suggests  they 
may  be  due  to  impurity  modes.  High-resolution  spectra  of  the  £!jlgh,  A[(LO),  and  ii^LO)  modes  were  recorded, 
and  an  analysis  of  the  anharmonicity  and  lifetimes  of  these  phonons  is  carried  out.  The  E2 gh  mode  displays  a 
visibly  asymmetric  line  shape.  This  can  be  attributed  to  anharmonic  interaction  with  transverse  and  longitudi¬ 
nal  acoustic  phonon  combinations  in  the  vicinity  of  the  K  point,  where  the  two-phonon  density  of  states 
displays  a  sharp  edge  around  the  £6lgh  frequency.  The  temperature  dependence  of  the  linewidth  and  frequency 
of  the  £6lgh  mode  is  well  described  by  a  perturbation-theory  renormalization  of  the  harmonic  E i)18*1  frequency 
resulting  from  the  interaction  with  the  acoustic  two-phonon  density  of  states.  In  contrast,  the  A^LO)  and 
£i(LO)  frequencies  lie  in  a  region  of  nearly  flat  two-phonon  density  of  states,  and  they  exhibit  a  nearly 
symmetric  Lorentzian  line  shape  with  a  temperature  dependence  that  is  well  accounted  for  by  a  dominating 
asymmetric  decay  channel. 
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I.  INTRODUCTION 

The  interest  in  ZnO  has  been  renewed  by  its  potential 
applications  in  transparent  electronics  and  UV  optoelectronic 
devices.  Its  large  exciton  binding  energy  (60  meV)  makes  it 
highly  attractive  for  applications  as  blue  and  ultraviolet  light 
emitters,  where  it  may  become  an  advantageous  alternative 
to  GaN.1  Good  knowledge  of  the  vibrational  properties  of 
this  material  is  essential  to  understand  transport  properties 
and  phonon  interaction  with  the  free  carriers,  both  of  which 
have  great  impact  on  optoelectronic  device  performance.  To¬ 
day  this  is  made  possible  by  the  availability  of  high-quality, 
large-area  bulk  ZnO  single  crystals.  In  polar  semiconductors, 
carriers  excited  high  in  the  conduction  band  relax  toward 
their  ground  state  mainly  by  Frohlich  interaction  with  the 
longitudinal  optical  phonons.  Thus,  the  dynamics  of  the  pho¬ 
non  population  strongly  affects  the  performance  of  high¬ 
speed  optoelectronic  devices.2  Raman  scattering  measure¬ 
ments  are  well  suited  to  obtain  information  about  the  sample 
quality  in  wide  band  gap  semiconductors3-5  as  well  as  to 
analyze  more  specific  aspects  of  lattice  dynamics,  such  as 
isotopic  effects6-7  and  phonon  lifetimes.8,9  The  temperature 
dependence  of  the  phonon  modes  has  also  been  used  to  de¬ 
termine  the  local  temperature  of  a  GaN  diode  as  a  function  of 
the  operating  voltage,10  and  therefore  Raman  scattering  has 
proven  to  be  well  suited  to  monitor  the  local  temperature 
during  device  operation. 


Wurtzite-type  ZnO  belongs  to  the  space  group  C^v  with 
two  formula  units  in  the  primitive  cell.  The  zone-center  op¬ 
tical  phonons  can  be  classified  according  to  the  following 
irreducible  representations:  F(lp[=A ,  +E]  +2E’2  +  2/i1.  The  Bl 
modes  are  silent  modes,  the  A1  and  Ei  modes  are  polar 
modes  and  are  both  Raman  and  infrared  active,  whereas  the 
E2  modes  are  nonpolar  and  Raman  active  only.  Polarized 
right-angle  Raman  scattering  of  ZnO  was  measured  by  Da- 
men  et  al.n  back  in  1965.  Arguello  et  al.n  reported  Raman 
spectra  of  ZnO  obtained  in  a  variety  of  geometries,  which 
allowed  them  to  identify  all  Raman  active  modes  from  their 
selection  rules.  A  few  years  later,  a  detailed  study  of  the 
resonant  behavior  of  the  high  frequency  E2  mode  ( /'’t,"gh) , 
A ,(TO),  and  £'1(TO,LO)  modes  as  well  as  several  second- 
order  features  was  presented  by  Calleja  and  Cardona.13  Data 
on  temperature  dependence  of  ZnO  phonons  also  date  back 
to  the  late  1970s.14  In  the  past  few  years,  the  availability  of 
bulk  ZnO  crystals  of  very  high  optical  quality  and  the  tech¬ 
nological  relevance  of  this  material  have  spurred  the  realiza¬ 
tion  of  fundamental  and  more  detailed  studies  on  ZnO  lattice 
dynamics.  However,  our  knowledge  of  the  lattice  dynamics 
of  ZnO  is  still  rather  limited,  as  few  experimental  data  are 
available  on  phonon  dispersion  in  this  material.  In  fact,  de¬ 
tailed  measurements  of  the  optical  phonon  branch  by  inelas¬ 
tic  neutron  scattering  are  still  lacking.  Only  recently  have  ab 
initio  density-functional  theory  (DFT)  calculations  of  lattice- 
dynamical  properties  been  published,15  which  have  made  a 


1098-0 12 1/2007/75(16)/!  65202(11) 


165202-1 


©2007  The  American  Physical  Society 


CUSCO  et  al. 


PHYSICAL  REVIEW  B  75,  165202  (2007) 


TABLE  I.  Symmetry  allowed  Raman  modes  for  the  scattering 
geometries  considered  in  this  work. 


Scattering 

geometry 

E2 

Raman  active  modes 

A^LO)  AJ  (TO)  Ei(LO) 

£j(TO) 

z(xx)z 

A 

A 

z(xy)z 

A 

x(yy)x 

A 

A 

x(zy)x 

A 

x(zy)y 

A 

A 

more  accurate  interpretation  of  the  experimental  data  pos¬ 
sible.  Thus,  the  unusual  variations  of  the  E^gh  phonon  line- 
width  with  the  isotopic  composition  and  with  pressure  have 
been  accounted  for  by  fine  features  of  the  phonon  density  of 
states  (DOS)  revealed  by  the  DFT  calculations.8  Recently, 
the  low-frequency  E2  phonon  lifetime  has  been  measured  in 
ZnO  by  means  of  impulsive  stimulated  Raman  scattering  ex¬ 
periments,  and  it  has  been  found  to  be  longer  than  200  ps  at 
low  temperature.16  It  has  been  suggested  that  such  long-lived 
optical  phonons  could  have  applications  in  quantum  cryptog¬ 
raphy  and  quantum  computing.16 

In  this  work  we  present  a  study  of  the  first-  and  second- 
order  Raman  spectra  of  ZnO  and  their  dependence  on  tem¬ 
perature.  Second-order  Raman  features  are  discussed  in  the 
light  of  recently  available  DFT  calculations  of  phonon  DOS. 
To  analyze  the  temperature  dependence  of  optical  phonons  in 
ZnO,  high  spectral  resolution  measurements  of  the  E 6lsh, 
A  |  (LO),  and  /T,  (LO)  modes  have  been  carried  out  at  different 
temperatures.  A  visible  line-shape  asymmetry  is  observed  for 
the  Zsi}18'1  mode,  which  can  be  explained  by  the  features  of  the 
two-phonon  DOS  around  the  E2lgh  frequency.  We  discuss  the 
temperature  effects  on  the  phonons  in  terms  of  third-order 
anharmonic  decay.  Whereas  the  behavior  of  the  A  |(LO)  and 
£[(LO)  modes  is  well  described  by  a  dominant  decay  chan¬ 
nel  localized  around  a  single  set  of  frequencies,  the  descrip¬ 
tion  of  the  E^'gh  temperature  dependence  requires  a  more 
detailed  analysis  that  takes  into  account  the  two-phonon 
DOS. 


of  the  spectrometer  with  100  /am  slits,  corresponding  to  a 
spectral  slit  width  of  ~2.5  cm-1.  The  triple  additive  configu¬ 
ration  with  50  /J,m  slit  was  utilized  to  obtain  high  resolution 
spectra  of  the  E2'gh,  A ,  ( LO) ,  and  Z8]  (LO)  modes.  The  spec¬ 
trometer  spectral  broadening  for  this  configuration  was 
checked  by  measuring  the  \  =  546.07  nm  emission  line  of 
Hg,  which  was  fitted  by  a  Gaussian  line  shape  with  a  full 
width  at  half  maximum  (FWFIM)  of  rc=0.44  cm-1.  Then, 
following  Ref.  8,  the  instrumental  broadening  was  estimated 
from  rL=r-r“/r,  where  T  and  TL  are,  respectively,  the 
FWHM  measured  on  the  Raman  spectrum  and  the  intrinsic 
FWHM  of  the  mode.  Using  this  relation,  we  find  that  the 
instrumental  broadening  correction  for  the  narrowest  mode 
studied  (the  E^gb  mode  at  80  K)  amounts  to  less  than  1%  of 
the  measured  FWHM,  and  therefore  we  can  assume  that  the 
high-resolution  Raman  spectra  reflect  the  intrinsic  FWHM  of 
the  modes. 

A  LN2  cryostat  was  used  to  vary  the  sample  temperature 
from  80  to  380  K.  Within  this  temperature  range,  the  spectra 
were  acquired  in  macrocamera  configuration.  To  extend  the 
temperature  range  up  to  750  K,  we  used  a  Linkam  high- 
temperature  stage  and  the  measurements  were  performed  us¬ 
ing  a  confocal  microscope  with  a  long  working-distance  X  80 
objective  protected  by  a  cooling  jacket.  Both  sets  of  mea¬ 
surements  yield  consistent  results  in  the  overlapping  tem¬ 
perature  range  between  290  and  380  K.  For  each  measure¬ 
ment,  the  temperature  was  stabilized  for  10  min  before 
acquiring  a  spectrum. 


III.  RESULTS  AND  DISCUSSION 
A.  First-  and  second-order  Raman  scattering  of  ZnO 

First,  we  present  an  overview  of  the  first-order  Raman 
scattering  in  ZnO,  and  compare  our  results  with  previously 
reported  data.  We  have  performed  polarized  Raman  scatter¬ 
ing  measurements  in  different  scattering  configurations  to 
identify  all  Raman  active  modes  of  ZnO  (see  Table  I).  The 
Raman  tensors  for  each  of  these  modes  are  given  by 
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II.  EXPERIMENT 

The  experiments  were  performed  on  a  high-quality  ZnO 
single  crystal  obtained  by  the  hydrothermal  growth  method. 
Growth  details  as  well  as  x-ray  and  photoluminescence  (PL) 
characterization  of  ZnO  crystals  obtained  by  this  method 
have  been  published  elsewhere.17  After  chemomechanical 
polish,  PL  spectra  obtained  at  2.1  K  revealed  emission  from 
excitonic  transitions  with  narrow  line  shapes,  indicating  the 
high  quality  of  these  hydrothermally  grown  samples.17 

The  Raman  spectra  were  excited  with  the  514.5  nm  line 
of  an  Ar+  laser  in  different  scattering  geometries  (see  Table 
I).  The  scattered  light  was  analyzed  by  means  of  a  Jobin- 
Yvon  T64000  triple  spectrometer  equipped  with  a 
LN2-cooled  charge-coupled  device  (CCD)  detector.  Full- 
range  spectra  were  obtained  in  the  subtractive  configuration 
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For  polar  modes,  the  coordinate  in  parentheses  denotes 
the  direction  of  phonon  polarization.  As  derived  from  the 
Raman  tensors  of  the  respective  modes,  all  Raman  active 
modes  can  be  observed  in  backscattering  geometry  except 
for  the  Zs^LO)  mode,  for  which  only  the  polarizability  com¬ 
ponents  involving  z  and  the  direction  of  phonon  polarization 
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FIG.  1 .  Room  temperature  first-order  Raman  spectra  of  ZnO  in 
scattering  geometries  corresponding  to  Ax+E2  and  Ex  symmetries. 

change  during  the  vibration.  The  symmetry-allowed  Raman 
modes  for  the  scattering  geometries  considered  in  this  work 
are  summarized  in  Table  I  for  convenience.  Although  the 
Zs^LO)  mode  was  reported  in  the  forbidden  configuration 
x(zz)x  under  406.7-nm  excitation  owing  to  intraband 
Frohlich  interaction,13  the  unambiguous  observation  of  the 
Zs^LO)  mode  requires  a  scattering  configuration  with  a  non¬ 
vanishing  x(zy)y  component. 

Figure  1  shows  polarized  first-order  Raman  spectra  of  the 
ZnO  sample  obtained  at  room  temperature  in  four  different 
scattering  configurations  that  allow  the  observation  of  all  the 
Raman  active  modes.  Our  results  are  consistent  with  previ¬ 
ous  Raman  studies. 11-13  The  z(xx)z  spectrum  ( E2 + A ,  sym¬ 
metry)  is  dominated  by  the  two  intense,  sharp  E2  modes  at  99 
and  438  cm-1.  The  low-frequency  E2  mode,  involving 
mainly  Zn  motion,  displays  an  extremely  narrow  linewidth, 
well  below  the  slit  width  of  the  subtractive  configuration  of 
the  spectrometer  even  for  the  highest  temperature  studied. 
The  E^'gh  mode  displays  a  clear  asymmetry  toward  low  fre¬ 
quencies  that  we  shall  discuss  in  more  detail  in  the  next 
section.  In  contrast  with  previous  works,  where  the  A ,  ( LO) 
mode  could  not  be  detected  for  excitation  wavelengths 
longer  than  406.7  nm,L'  we  can  detect  a  weak  ,4 ,  ( LO )  mode 
at  574  cm-1.  The  prominent  feature  at  333  cm-1  corresponds 
to  second-order  scattering  and  will  be  discussed  below.  In  the 
x(yy)x  spectrum  (E2+A ,  symmetry),  in  addition  to  the  strong 
E2  modes  a  new  peak  appears  at  378  cm-1,  which  can  be 
assigned  to  the  A^TO)  mode.  An  additional  peak  emerges  at 
410  cm-1  in  the  x(zy)x  spectrum  (Et  symmetry),  which  is 
assigned  to  the  Zs^TO)  mode.  Finally,  the  E ,  ( LO)  mode  can 
be  observed  at  590  cm-1  in  the  x(zy)y  spectrum.  Therefore, 
the  highest-frequency  first-order  mode  of  ZnO  corresponds 
to  the  Z),  (LO)  mode,  which  at  zone  center  is  found  at  a 
frequency  16  cm-1  higher  than  that  of  the  A^LO)  mode.  It  is 
worth  noting  that  DFT  calculations  yield  inaccurate  results 
regarding  the  Ti,  (LO)  and  A  |  (LO)  modes,  as  they  give  a  zone 
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center  Zq(LO)  frequency  4  cm-1  lower  than  that  of  the 
A ,  ( LO)  mode.15  This  fact  should  be  taken  into  account  when 
discussing  the  Raman  features  of  ZnO  in  terms  of  the  calcu¬ 
lated  phonon  DOS.  The  observed  first-order  Raman  frequen¬ 
cies  are  reported  in  Table  II. 

A  residual  intensity  of  the  forbidden  E2  modes  can  be 
seen  in  the  x(zy)x  and  x(zy)y  spectra  that  could  arise  from  a 
slight  misalignment  of  the  sample  introducing  a  small  x  com¬ 
ponent  in  the  z  polarization  direction  or  from  a  lower  quality 
of  the  sample  polishing  in  the  m  prismatic  face  resulting  in  a 
relaxation  of  the  selection  rules. 

Whereas  wave-vector  conservation  restricts  the  phonons 
involved  in  first-order  Raman  scattering  to  those  with  k  ~  0, 
phonons  from  the  entire  Brillouin  zone  take  part  in  second- 
order  Raman  scattering.  Therefore,  second-order  spectra  usu¬ 
ally  display  feature-rich  structures,  which  are  determined,  on 
the  one  hand,  by  the  phonon  DOS  and,  on  the  other  hand,  by 
the  selection  rules  of  the  two-phonon  scattering  processes. 
Selection  rules  for  two-phonon  Raman  scattering  in  crystals 
with  the  wurtzite  structure  were  reported  by  Siegle  et  al. 18  As 
a  general  rule,  it  was  found  that  overtones  always  contain  the 
representation  A ,  whereas  combinations  of  phonons  belong¬ 
ing  to  different  representations  do  not  contain  the  A:  repre¬ 
sentation.  According  to  DFT  calculations15  the  phonon  DOS 
of  ZnO  presents  a  frequency  gap  between  acoustic  and  opti¬ 
cal  modes  that  extends  from  270  to  410  cm-1.  The  second- 
order  spectra  may  then  be  divided  into  three  regions:  (i)  the 
low-frequency  region  (approximately  160-540  cm-1)  domi¬ 
nated  by  acoustic  overtones,  (ii)  the  high-frequency  region 
(820-1120  cm-1)  formed  by  optical  overtones  and  combina¬ 
tions,  and  (iii)  the  intermediate-frequency  region 
(540-820  cm-1)  where  optical  and  acoustic  phonon  combi¬ 
nations  occur. 

In  the  present  analysis,  we  use  the  phonon  dispersion  cal¬ 
culated  by  DFT  (Ref.  15)  and  the  symmetry  selection  rules 
for  two-phonon  Raman  scattering  reported  by  Siegle  et  al. 13 
to  aid  the  identification  of  the  second-order  features  of  the 
ZnO  Raman  spectrum.  Figure  2  shows  the  z(xx)z  polarized 
second-order  Raman  spectra  of  ZnO  ( A{+E2  symmetry)  for 
different  temperatures  in  the  80-300  K  range.  The  second- 
order  features  are  labeled  with  their  respective  frequencies 
on  the  RT  spectrum.  The  most  prominent  second-order  fea¬ 
tures  occur  in  the  high-frequency  region  and  correspond  to 
LO  overtones  and  combinations  involving  LO  modes.  The 
broad,  intense  peak  at  1158  cm-1,  which  is  found  between 
the  doubled  frequencies  measured  for  the  A ,  (LO)  and 
E]  ( LO)  modes,  contains  contributions  of  24 ,  ( LO)  and 
2/T|  (LO)  modes  at  the  I  point  of  the  Brillouin  zone,  and 
possibly  also  of  2LO  scattering  by  mixed  modes  from  the 
rather  flat  bands  along  the  A-L-M  line.  The  weaker  peak  at 
1 105  cm-1  can  be  attributed  to  2LO  at  ZZ  and  K.  Note  that  the 
DFT  calculations,  in  addition  to  giving  the  wrong  ordering 
for  the  E ,  ( LO )  and  4 ,  ( LO )  modes,  tend  to  underestimate  the 
LO  frequencies.  Thus,  the  Raman  peaks  in  this  region  occur 
at  higher  frequencies  than  predicted  by  DFT  results.  As  men¬ 
tioned  above,  all  these  overtones  contain  the  Al  representa¬ 
tion,  and  furthermore,  they  occur  in  a  region  with  high  two- 
phonon  DOS.  Two  weak  shoulders  can  be  observed  in  the  RT 
spectrum  at  1072  and  1044  cm-1,  which  become  more  vis- 
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TABLE  II.  Room  temperature  frequencies  and  symmetries  of  the  first-  and  second-order  Raman  spectra 
observed  in  ZnO  and  their  assignments.  Our  results  are  compared  with  previous  data  in  Ref.  13.  Parentheses 
indicate  symmetries  that  although  being  present  in  the  spectra  display  a  much  lower  intensity  than  the 
dominant  one. 


Frequency  (i 

Our  data 

cm  :) 

Ref.  13 

Symmetry 

Process 

Brillouin  zone 
points/lines 

99 

101 

e2 

Tp  lOW 

1 

203 

208 

Au(E2) 

2TA;  2f?2°w 

L,M,H\ r 

284 

^1 

jghigh_glow 

1 

333 

332 

/-high  plow 
^ 2  -c2 

1 

378 

380 

^1 

A!  (TO) 

1 

410 

408 

E\ 

Ei  (TO) 

1 

438 

437 

e2 

£>high 

1 

483 

^1 

2LA 

M-K 

536 

541 

^1 

2Blow;  2LA 

T \L,M,H 

574 

574 

^1 

A  \  (LO) 

1 

590 

584 

Ei 

Ei(LO) 

1 

618 

^1 

TA+TO 

H,M 

657 

Ei,E2 

TA+LO 

L,H 

666 

^1 

TA+LO 

M 

700 

^1 

LA+TO 

M 

723 

^1 

LA+TO 

L-M 

745 

^1 

LA+TO 

L-M 

773 

^1 

LA+TO 

M  ,K 

812 

^1 

LA+LO 

L,M 

980 

990 

^1 

2TO 

L-M-K-H 

1044 

^1 

TO+LO 

AM 

1072 

1080 

^1 

TO+LO 

M  ,L 

1105 

^1 

2LO 

H,K 

1158 

1160 

^1 

2A1(LO),2£1(LO);  2LO 

t-a-l-m 

ible  at  low  temperature.  We  assign  the  mode  at  1072  cm-1  to 
TO+LO  combinations  at  the  M  and  L  points,  and  the  mode 
at  1044  cm-1  to  TO+LO  combinations  at  the  A  and  H  points. 
A  weak  broad  band  can  be  observed  at  about  980  cm-1, 
which  can  be  attributed  to  TO  overtones.  It  is  noteworthy 
that  the  scattering  efficiency  of  processes  involving  LO 
modes  is  much  higher  than  those  involving  TO  modes.  As 
reported  by  Calleja  and  Cardona,13  the  2LO  scattering  is 
resonantly  enhanced  already  at  photon  energies  ~1  eV  be¬ 
low  the  fundamental  gap. 

A  multiplicity  of  structures  appear  in  the  low-  and 
intermediate-frequency  regions  of  the  spectra,  where  acous¬ 
tic  overtones  and  acoustic-optical  combinations  occur.  To  aid 
their  identification  by  analyzing  their  symmetry,  we  have  re¬ 
corded  polarized  Raman  spectra  in  several  scattering  geom¬ 
etries.  In  Fig.  3  we  display  the  At+E2,  A1,  E2,  and  El  com¬ 
ponents  of  the  second-order  spectrum.  The  most  intense  peak 
in  the  intermediate-low-frequency  region  is  observed  at 
333  cm-1.  This  mode  had  been  previously  assigned  to  trans¬ 
verse  acoustic  overtone  scattering  at  M. 13,19  However,  the 
temperature  dependence  of  its  Raman  intensity  (see  Fig.  2) 
clearly  indicates  that  this  is  a  difference  mode.  The  fre¬ 


quency  of  this  mode  is  in  good  agreement  with  the  difference 
between  the  E^lsh  and  E2m'  frequencies  measured  in  our 
sample.  Using  the  notation  of  Ref.  18,  the  E^'^—E1™  differ¬ 
ence  contains  symmetries  Tg^rgDA!,^-  This  is  consistent 
with  the  polarized  measurements  shown  in  Fig.  3,  which 
indicate  that  the  symmetry  of  the  333  cm-1  mode  is  predomi¬ 
nantly  A  ] ,  with  a  smaller  E2  component  and  an  even  smaller 
El  component.  According  to  the  calculated  phonon  disper¬ 
sion  relations,15  the  333-cirr1  mode  could  also  contain  con¬ 
tributions  from  [TO-TA]A£W  differences.  These  differ¬ 
ences,  however,  yield  higher  frequencies  and  all  of  them 
contain  the  El  symmetry.  To  demonstrate  the  origin  of  the 
333-cm_1  mode,  we  plot  in  Fig.  4  the  Raman  intensity  of  this 
mode  versus  temperature,  normalized  to  its  intensity  at  room 
temperature,  and  we  compare  it  with  the  statistical  occupa¬ 
tion  factor  for  the  E2lgh - E2"  difference  mode  (dashed  line), 
which  is  given  by20 

p(T)  =  [l+n(Ei?sh)]n(E'™),  (2) 

where  n(E)  =  \exp(E/ kBT)~  1]_1  is  the  Bose-Einstein  distri¬ 
bution  function.  For  comparison,  we  also  plot  (dotted  line) 
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FIG.  2.  Second-order  Raman  spectra  of  ZnO  obtained  in  z(xx)z 
configuration  for  temperatures  in  the  80-300  K  range. 

the  statistical  occupation  factor  for  overtones  of  TA  phonons 
with  ftcoXA~  160  cm”1  around  the  M  point,  p(T)  =  \  I 
+n(ftcoXA)]2,  as  suggested  in  Refs.  13  and  19.  Clearly,  the 
statistical  occupation  factor  for  the  difference  mode  [Eq.  (2)] 
accounts  for  the  observed  temperature  dependence  of  the  in¬ 
tensity,  whereas  the  corresponding  factor  for  the  overtones 
deviates  significantly  from  the  experimental  points. 

Another  prominent  peak  is  observed  at  203  cm-1,  which 
as  seen  from  Fig.  3  exhibits  Ay  symmetry  and  thus  can  be 
attributed  to  a  TA  overtone.  In  Ref.  15  this  mode  was  as¬ 
signed  to  2TAl,  and  in  Ref.  13  to  2 E1™  with  possible  con¬ 
tributions  of  2TA  at  the  M  point.  We  note  that  the  corre¬ 
sponding  phonon  branch  is  rather  flat  at  //,  so  we  speculate 
that  the  203  cm-1  peak  could  also  have  a  significant  contri¬ 
bution  of  2TA  scattering  at  the  H  point. 

A  weak  peak  can  be  seen  at  284  cm-1  in  the  room  tem¬ 
perature  spectrum,  which  disappears  at  lower  temperatures 
(see  Fig.  2),  and  therefore  we  assign  it  to  a  difference  mode. 
As  can  be  seen  from  Fig.  3,  this  mode  has  Ay  symmetry. 
Several  mode  differences  yield  frequencies  close  to  the  ob¬ 
served  peak:  A1(LO)-,B11ow,  £1(FO)-B11ow,  B\ieh-Bl™,  and 
[LO-LAj/.  However,  some  of  them  have  symmetries  which 
are  incompatible  with  the  spectra  of  Fig.  3.  The  symmetry  of 
these  differences  are,  respectively,  I  }  ®  r3 D  T5 
®  T3 IDE2,  r3®r3DA1,  and  Ly®  LyDAy  ,Ey,E2.1&  Therefore, 
the  first  two  possibilities  can  be  excluded,  and  the  absence  of 
Ey  and  E2  component  in  this  peak  suggests  that  it  can  be 
primarily  assigned  to  Z?ljllgh-5110w.  The  very  weak  structure  de¬ 
tected  at  483  cm-1,  which  as  seen  from  Fig.  3  exhibits  Ay 
symmetry,  can  be  assigned  to  LA  overtones  along  M-K, 
where  phonon  dispersion  is  nearly  flat  and  hence  the  DOS  is 
very  high.  The  peak  at  536  cm-1  is  clearly  of  Ay  symmetry 
and  can  be  attributed  to  2/illow  and  LA  overtones  along  L-M 
and  H. 

Above  this  frequency,  we  find  acoustic  and  optical  com¬ 
binations.  The  peak  at  618  cm-1  could  be  tentatively  as¬ 
signed  to  TA+TO  at  the  H  and  M  points.  TA+LO  combi¬ 
nations  give  rise  to  the  doublet  observed  at  657-666  cm-1. 
The  lower  component  of  the  doublet  has  E2  and  Ey  symme¬ 
try,  whereas  the  symmetry  of  the  higher  component  is  basi¬ 


FIG.  3.  Second-order  Raman  spectra  of  ZnO  obtained  for  vari¬ 
ous  scattering  geometries.  From  top  to  bottom,  z(xx)z,  x(zz)x, 
z(xy)z,  and  x(zy)x. 

cally  Ay.  The  highest  DOS  for  combinations  in  this  fre¬ 
quency  range  is  found  at  the  L,  M,  and  H  points.  Both  at  the 
L  and  H  points  the  TA+LO  combination  has  symmetry 
Ey,E2,  whereas  at  the  M  point  combinations  of  the  type  M, 
® M,DAy,E2  occur.  This  suggests  the  assignment  of  the 
657  cm-1  mode  to  TA+LO  combinations  at  the  L  and  H 
points  and  the  666  cm"1  mode  to  TA+LO  combinations  at 
the  M  point.  Low  intensity  modes  of  mainly  Ay  symmetry 
are  detected  at  700,  723,  745,  773,  and  812  cm-1,  which  can 
be  attributed  to  acoustic  and  optical  combinations.  Taking 
into  account  the  predominantly  Ay  symmetry  of  all  these 
modes,  we  tentatively  assign  the  peak  at  700  cm-1  to  LA 
+TO  combinations  at  the  M  point,  and  the  peak  at  812  cm"1 
to  LO  +  LA  combinations  at  the  L  and  M  points.  The  frequen¬ 
cies  of  the  second-order  Raman  features  observed  at  room 
temperature  are  collected  in  Table  II. 

It  is  noteworthy  that,  in  addition  to  the  second-order  peaks 
discussed  above,  several  sharp  peaks  appear  and  grow  in 
intensity  as  the  temperature  decreases.  These  peaks  are 


T  (K) 


FIG.  4.  Intensity  of  the  second-order  peak  observed  at 
=  333  cm"1  as  a  function  of  temperature.  Dashed  line,  occupation 
number  factor  for  the  difference  mode;  dotted  line,  idem 

for  TA  overtones  around  the  M  point. 
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marked  by  arrows  in  the  80  K  spectrum  shown  in  Fig.  2.  The 
weak  peak  at  380  cm-1  appears  at  the  same  frequency  as  the 
A  [(TO)  mode  in  a  x(yy)x  spectrum  recorded  at  80  K,  so  we 
assign  it  to  a  leaky  A  [(TO)  mode.  While  the  peak  at 
589  cm-1  is  close  to  the  /;  |  (LO)  mode,  its  frequency  is  about 
4  cm-1  lower  than  that  measured  for  the  £[(LO)  mode  at 
80  K  in  the  allowed  x(zy)y  scattering  configuration.  This 
peak  could  be  related  to  local  vibrational  modes  associated 
with  intrinsic  lattice  defects,  similarly  to  the  peak  that  was 
observed  at  582  cm-1  in  doped  GaN, 3,19,21  which  was  con¬ 
jectured  to  be  associated  with  intrinsic  lattice  defects  favored 
or  activated  by  doping.21  The  intensity  enhancement  and 
sharpening  with  decreasing  temperature  observed  for  this 
peak  as  well  as  for  that  at  142  cm-1  clearly  suggest  a  local 
mode  origin  for  these  peaks.  Although  Al,  Fe,  Si,  Li,  and  K 
impurities  were  detected  by  glow  discharge  mass  spectrom¬ 
etry  in  hydrothermal  ZnO  samples,  their  concentration  was 
found  to  be  in  the  few  ppm  range,17  well  below  the  detection 
limit  of  Raman  spectroscopy.  Thus,  we  speculate  that  these 
peaks  are  related  to  intrinsic  lattice  defects,  which  may  be 
favored  by  the  growth  conditions.  These  crystals  present  a 
strong  visible  luminescence  band,  which  supports  the  exis¬ 
tence  of  intrinsic  point  defects.  On  the  other  hand,  surface 
effects  can  be  ruled  out  because  of  the  deep  penetration  of 
the  excitation  light.  Further  work  should  be  carried  out  to 
determine  the  origin  of  these  Raman  peaks. 


B.  Anharmonic  phonon  decay  and  phonon  lifetimes 
1.  The  E\ )lgh  mode:  Fermi  resonance  model 

In  the  discussion  of  the  first-order  Raman  spectra  of  ZnO 
presented  in  the  preceding  section,  we  noted  that  the  E^lsh 
mode  exhibits  a  visibly  asymmetric  line  shape  with  a  low- 
frequency  tail.  This  is  quite  apparent  from  Fig.  1  when  one 
compares  the  line  shapes  of  the  ,E2lgh  and  A[(TO)  modes. 
Given  that  the  A  [(TO)  line  shape  is  fairly  symmetric,  the 
asymmetry  of  the  £1)1Eh  mode  cannot  be  ascribed  to  lattice 
disorder.  Furthermore,  isotopic  broadening  is  negligible  for 
the  Zi2lgh  mode  since  it  mainly  involves  O  motion,  and  O  is 
nearly  isotopically  pure.  The  line-shape  broadening  is  then 
mostly  determined  by  anharmonic  phonon-phonon  interac¬ 
tions.  These  can  result  in  strongly  distorted  peaks  when  reso¬ 
nant  interaction  with  a  band  of  second-order  combinations 
takes  place  (Fermi  resonance),  as  is  the  case,  for  instance,  for 
the  GaP  TO  mode,22  23  where  the  presence  of  van  Hove-type 
singularities  in  the  DOS  of  the  TA+LA  combination  band 
gives  rise  to  a  highly  asymmetric  TO  mode,  which  develops 
a  side  band  at  high  pressure.22  A  similar  situation  occurs  for 
the  £i)Ish  mode  of  ZnO,  as  its  frequency  lies  close  to  a  ridge¬ 
like  structure  of  the  two-phonon  DOS  corresponding  to  TA 
+LA  combinations  in  the  vicinity  of  the  K  point.8 

To  evaluate  this  effect  for  the  L,'sh  phonon  of  ZnO  we 
consider  the  phonon  self-energy  2(w)  =  A(w)-/Y(co),  which 
expresses  the  renormalization  of  the  bare  harmonic  fre¬ 
quency  of  the  mode,  oj0,  as  a  result  of  the  interaction  with 
other  phonons.  To  the  second  order  in  perturbation  theory, 
the  cubic  anharmonicity  effect  on  the  imaginary  part  of  the 
self-energy  is  given  by24 


1  877 

YM  =  —T  2  I  V3(qj'i,-  q/2)|2[(«i  +  n 2  +  1)<S(«1  +  0)2  -  co) 

h~  1/172 

+  (n2-nl)S(wl  -  co2-  to)],  (3) 

where  ,-<1/2)  is  the  third-order  coefficient  in  the  ex¬ 

pansion  of  the  lattice  potential  in  normal  coordinates  and  nx 
and  n2  are  the  Bose-Einstein  occupation  factors  for  the 
modes  <w(q,y'[)  and  <w(-qj2),  respectively  (q  is  the  phonon 
wave  vector  and  jt  labels  the  phonon  branch).  Assuming 
slowly  varying  Bose-Einstein  factors  and  U3  coefficients  in 
the  small  energy  range  considered,  the  imaginary  part  of  the 
self-energy  is  proportional  to  the  two-phonon  DOS  and  can 
be  written  as 

Y(w)  =  |  V3 12(  1  +  77  [  +  n2)p+(w)  +  |Yj|  2(n2-U[)p"(w), 

(4) 

with  p+(co)  and  p“(w)  the  two-phonon  sum  and  difference 
DOS,  and  Uj  and  L]  the  effective  third-order  anharmonic 
constants  for  decay  into  phonon  sums  and  phonon  differ¬ 
ences,  respectively.  Since  real  and  imaginary  parts  of  the 
phonon  self-energy  are  related  by  a  Kramers-Kronig  trans¬ 
formation,  the  corresponding  real  part  of  the  self-energy 
A  (&>)  can  be  evaluated  as 

2  r  o)' 

A(to)  =  -  ~V  2Y((o')dco' .  (5) 

77  JQ  to  -w“ 

This  yields  a  frequency-dependent  renormalization  of  the 
phonon  energy,  which  has  to  be  added  to  the  constant  shift  in 
frequency  due  to  the  thermal  expansion  of  the  lattice  A0. 
This  latter  contribution  can  be  written  as25 

A0  =  -co0y[  [, ac(T')  +  2aa(T')-\dT \  (6) 

Jo 


where  ac  and  aa  are  the  linear  thermal  expansion  coefficients 
along  directions  parallel  and  perpendicular  to  the  c  axis,  re¬ 
spectively,  and  y  is  the  Griineisen  parameter  of  the  /sijlgh 
mode,  which  we  take  as  y=  2.02. 15  We  have  derived  the  ther¬ 
mal  expansion  coefficients  for  ZnO  from  the  temperature  de¬ 
pendence  of  the  lattice  parameters  determined  by  x-ray  dif¬ 
fraction  by  Iwanaga  et  al.26  Although  these  data  are  only 
available  for  7’>  300  K,  we  have  extrapolated  the  results 
down  to  T=  200  K.  Below  this  temperature,  thermal  expan¬ 
sion  usually  falls  rapidly  and  becomes  negligible  at  T 
~  50  K,27  and  consequently  the  extrapolation  of  high- 
temperature  values  is  not  a  good  approximation.  Thus,  as  a 
more  accurate  approximation  for  the  thermal  expansion  co¬ 
efficients  for  T <  200  K  we  have  used  the  calculated  varia¬ 
tion  for  ZnSe,27  scaled  to  fit  smoothly  with  the  curve  derived 
from  experimental  results  for  ZnO.  For  small  self-energy 
corrections  compared  to  the  Raman  frequency,  the  line  shape 
of  the  Raman  peak  is 


/(to)  °c 


_ YM _ 

[too  +  Aq  +  A(to)  —  to]-  +  Y(to)- 


(7) 


If  the  frequency  dependence  of  the  phonon  self-energy  is 
neglected,  the  Raman  line  shape  given  by  Eq.  (7)  reduces  to 


165202-6 


TEMPERATURE  DEPENDENCE  OF  RAMAN  SCATTERING... 


PHYSICAL  REVIEW  B  75,  165202  (2007) 


a  Lorentzian  function.  In  that  case,  the  imaginary  part  of  the 
phonon  self-energy  gives  directly  the  half  width  of  the  cor¬ 
responding  line  shape  Y  =  |r.  Note,  however,  that  because  of 
the  frequency  dependence  of  the  self-energy,  the  line  shape 
given  by  Eq.  (7)  may  deviate  substantially  from  a  Lorentzian 
line  shape  and  the  comparison  between  the  Raman  linewidth 
and  the  theoretical  lifetime  requires  a  detailed  analysis  that 
involves  a  calculation  of  both  the  real  and  the  imaginary  part 
of  the  self-energy.28  We  use  Eq.  (7)  with  the  frequency- 
dependent  self-energy  in  order  to  simulate  the  Raman  spec¬ 
trum  of  the  /it,llgh  mode  and  explain  the  observed  asymmetric 
line  shape. 

To  evaluate  the  imaginary  part  of  the  self-energy  we  have 
used  the  two-phonon  sum  DOS  reported  by  Serrano  et  al.,& 
shifted  by  -10  cm-1  as  explained  in  Ref.  15.  Following  the 
approach  taken  in  Refs.  6  and  23,  we  smooth  the  calculated 
DOS  in  order  to  avoid  spurious  artifacts  in  the  calculated  line 
shapes.  On  the  other  hand,  the  two-phonon  difference  DOS 
is  nearly  flat  around  the  £2lsh  frequency,  and  can  be  approxi¬ 
mated  by  a  constant  value  p~{u>)~  0.25  states/cm'1,  as  given 
in  Ref.  8.  According  to  the  calculated  phonon  dispersion  of 
ZnO,  the  most  likely  decay  channels  for  the  Z?!jlgh  mode  are 
the  decay  into  two  acoustic  modes  at  oq—^Ocm-1  and 
oj2~  250  cm-1,  and  into  a  difference  between  a  longitudinal 
optical  mode  at  o>2~5 50  cm-1  and  an  acoustic  mode  at  w, 
—  110  cm-1.8  Using  the  corresponding  Bose-Einstein  occu¬ 
pation  factors  in  Eq.  (4),  the  line  shape  of  the  /itllgh  mode  can 
be  calculated  from  Eqs.  (4)-(7)  with  two  adjustable  param¬ 
eters  |Vjj2  and  |Vj|2.  We  have  taken  <w0  =  440  cm-1  as  given 
by  the  ab  initio  calculations  reported  in  Ref.  15.  A  rigid  shift 
of  the  two-phonon  sum  DOS  has  been  included  to  take  into 
account  the  thermal  expansion  effects.  This  has  been  ob¬ 
tained  from  Eq.  (6)  with  the  average  Griineisen  parameters  of 
the  acoustic  modes  at  ftij  ~  190  cm-1  and  oj2~  250  cirT1, 
which,  based  on  the  values  reported  in  Ref.  15,  we  take  as 
y~  1.1. 

In  Fig.  5  we  plot  the  real  and  imaginary  parts  of  the  pho¬ 
non  self-energy  and  the  theoretical  line  shape  (solid  line)  for 
room  temperature  and  |Vjj2  =  39  cm-2  and  |Vjj2  =  5.8  cm-2. 
We  shall  discuss  the  determination  of  the  anharmonic  cou¬ 
pling  constants  below.  The  calculated  line  shape  has  been 
upshifted  by  2.3  cm-1  to  bring  it  into  agreement  with  the 
experimental  spectrum.  Such  frequency  shift  is  necessary  be¬ 
cause  of  the  low  bare  harmonic  frequency  given  by  the  ab 
initio  calculations,  which  was  estimated  to  be  about  0.7%  too 
low  in  Ref.  15.  Very  good  agreement  is  obtained  with  the 
high-resolution  Raman  spectrum  of  the  /-’o'gh  mode  (open 
circles).  As  can  be  seen  from  Fig.  5,  the  /T1) 1  s h  line  shape 
exhibits  an  asymmetric  broadening  toward  low  frequencies, 
where  it  deviates  significantly  from  the  Lorentzian  line  shape 
fitted  to  the  upper  frequency  half  of  the  peak  (dotted  line). 
The  asymmetry  in  the  E^lgh  line  shape  stems  from  the  strong 
frequency  dependence  of  the  imaginary  part  of  the  phonon 
self-energy,  which  is  plotted  in  Fig.  5  (dashed-dotted  line) 
against  the  left  axis  over  the  frequency  region  where  the  £2lgh 
mode  occurs.  Such  a  strong  variation  of  the  imaginary  part  of 
the  phonon  self-energy  in  the  vicinity  of  the  £2lgh  frequency 
has  been  previously  shown  to  give  rise  to  an  anomalous  be¬ 
havior  of  the  £!jI8h  linewidth  in  ZnO  crystals  with  different 


Frequency  (cm  1) 

FIG.  5.  High-resolution  spectrum  of  the  E2' 8,1  mode  of  ZnO  at 
room  temperature  (open  circles).  The  solid  line  is  the  theoretical 
line  shape  calculated  using  the  anharmonic  decay  model  described 
in  the  text.  For  comparison,  a  symmetric  Lorentzian  line  shape  fit¬ 
ted  to  the  high-frequency  side  of  the  peak  is  also  displayed  (dotted 
line).  The  real  and  imaginary  parts  of  the  phonon  self-energy  used 
in  the  calculations  are  plotted  as  dashed  and  dot-dashed  lines, 
respectively. 

isotopical  compositions  and  to  an  unusual  pressure  depen¬ 
dence  of  the  £^Igh  linewidth  vs  hydrostatic  pressure.8  In  the 
present  work,  we  show  that  this  self-energy  features  are  also 
responsible  for  the  line-shape  asymmetry  of  the  £!)lgh  Raman 
peak. 

The  determination  of  the  anharmonic  coupling  constants 
by  fitting  Eq.  (7)  to  a  single  Raman  spectrum  has  some  de¬ 
gree  of  arbitrariness,  as  many  combinations  of  |Vjj2  and  |Vj|2 
values  yield  similarly  good  visual  agreement  with  the  data. 
Instead  of  relying  on  a  single  Raman  spectrum  to  derive  the 
values  of  the  anharmonic  constants,  we  have  measured  the 
E^gh  linewidth  for  the  full  set  of  Raman  spectra  obtained  at 
different  temperatures  and  we  have  fitted  the  model  to  the 
temperature  dependence  of  the  E^gh  linewidth.  To  this  end, 
the  FWHM  obtained  from  the  theoretical  line  shape  [Eq.  (7)] 
rLS  was  tabulated  for  a  grid  of  |V^|2  and  |Vj|2  values.  The 
tabulated  values  were  used  to  interpolate  rLS(Y;|^|2,|Vj|2), 
which  was  fitted  by  the  Levenberg-Marquardt  method  to  the 
experimental  linewidths.  Including  a  constant  background 
contribution  to  the  linewidth  (T0)  to  account  for  defect  and 
impurity  scattering  improves  the  fit  significantly.  The  best  fit 
to  our  data,  which  is  displayed  in  Fig.  6,  is  obtained  for 
|V^[2  =  39  cm-2,  |Vj|2  =  3.8  cm-2,  and  r0=1.3  cm-1.  These  re¬ 
sults  confirm  that  the  decay  rate  for  phonon-difference  decay 
channels  is  considerably  lower  than  for  phonon-sum  decay 
channels,  and  thus  it  can  be  safely  neglected  at  low  tempera¬ 
tures. 

Estimates  of  the  anharmonic  squared  matrix  elements  for 
the  Efgh  mode  of  ZnO  were  obtained  in  Refs.  8  and  15  by 
using  the  expressions  for  the  imaginary  part  of  the  self¬ 
energy  [Eq.  (4)]  evaluated  at  the  frequency  of  the  mode  and 
adjusting  the  anharmonic  coefficients  to  reproduce  their  mea¬ 
sured  FWHM.  As  discussed  above,  when  the  frequency  de¬ 
pendence  is  neglected,  the  imaginary  part  of  the  self-energy 
corresponds  to  the  half  width  of  the  Lorentzian  line  shape, 
and  therefore  this  procedure  yields  values  for  the  anharmonic 
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FIG.  6.  Temperature  dependence  of  the  linewidth  (squares,  left 
axis)  and  frequency  (triangles,  right  axis)  of  the  E i)1®11  mode.  The 
dashed  line  is  a  fit  of  the  model  discussed  in  the  text  to  the  line- 
width  data.  The  dot-dashed  line  corresponds  to  the  temperature  de¬ 
pendence  of  the  frequency  given  by  the  model  for  the  same  set  of 
parameters. 


coefficients  which  are  actually  two  times  the  values  of  the 
corresponding  anharmonic  squared  matrix  elements  [|V^p 
and  |  V^|2  of  Eq.  (4)].  After  correcting  for  this  factor  of  2,  our 
values  for  the  anharmonic  squared  matrix  elements  are  some¬ 
what  higher  than  those  given  in  Refs.  8  and  15  (39  vs 
28  cm-1  for  |V^p  and  3.8  vs  3.1  cm-1  for  |V^p).  In  those 
works,  no  defect  and/or  impurity  contribution  to  the  line- 
width  was  considered.  If  we  do  not  include  the  constant  con¬ 
tribution  T0  in  our  fits,  the  resulting  |  V^\2  value  increases 
appreciably  to  compensate  for  the  defect  and/or  impurity 
broadening,  while  |V^p  turns  out  to  be  roughly  the  same. 
This  explains  the  value  |V^p=5.8  cm-2  used  to  obtain  the 
visual  fit  displayed  in  Fig.  5,  but  cannot  account  for  the 
lower  |V^P  value  reported  in  Ref.  8.  The  discrepancies  in  the 
|V^p  and  |V^p  values  mentioned  above  could  arise  from  the 
simplified  model  used  in  Refs.  8  and  15,  which  does  not  take 
into  account  the  actual  frequency  dependence  of  the  phonon 
self-energy.  Also,  one  should  bear  in  mind  that  the  higher- 
order  terms  neglected  in  Eq.  (3)  may  lead  to  deviations  in  the 
calculated  line  shapes  that  are  of  increasing  importance  for 
higher  temperatures.  This  could  explain  that  the  measured 
Fp'21'  frequency  falls  slightly  below  the  one  calculated  from 
the  model  for  7’>500  K  using  the  parameters  determined 
above  (see  Fig.  6). 29 

Figure  7  displays  high-resolution  E2lsh  spectra  recorded  at 
temperatures  in  the  80-750  K  range.  As  temperature  in¬ 
creases,  the  £pgh  peak  is  reduced  in  intensity  and  broadens 
significantly,  maintaining  a  visible  asymmetry  toward  lower 
frequencies.  Such  behavior  is  well  reproduced  by  the  line 
shapes  calculated  according  to  Eq.  (7),  which  are  plotted  in 
the  inset  of  Fig.  7.  Note  that  Eq.  (7)  does  not  contain  the 
defect  and/or  impurity  broadening  contribution  and  conse¬ 
quently  the  calculated  line  shapes  are  narrower  than  the  cor¬ 
responding  Raman  spectra,  which  are  further  broadened  by 
defects  and  impurities.  This  is  most  apparent  for  the  narrow¬ 
est  Raman  line  shape  measured  at  80  K. 


Frequency  (cm  p 


FIG.  7.  Raman  spectra  of  the  Epgh  mode  for  temperatures  in  the 
80-750  K  range.  Inset,  corresponding  line  shapes  calculated  using 
the  Fermi  resonance  model  described  in  the  text. 

2.  The  A  j(LO)  and  Ej(LO)  modes:  Generalized  Ridley  channel 

model 

The  A  |  (LO)  and  E](LO)  frequencies  of  ZnO  are  very 
close  to  each  other.  Unlike  the  Epsh  mode,  these  two  modes 
occur  at  frequencies  that  lie  in  a  plateau  of  the  two-phonon 
DOS  with  a  relatively  high  density  of  states.  Therefore  their 
line  shape  is  not  affected  by  variations  of  the  phonon  self¬ 
energy  and  a  simpler  model  can  be  used  to  analyze  their 
temperature  dependence.  Inspection  of  the  phonon  dispersion 
curves15  suggests  that  the  main  decay  channel  for  these 
modes  is  the  so-called  generalized  Ridley  channel.30  31  This 
involves  the  decay  of  the  longitudinal  optical  mode  into  a 
mode  of  the  transverse  acoustic  branches  and  a  mode  of  the 
transverse  optical  branches.  The  anisotropy  of  the  force  con¬ 
stants  in  the  wurtzite  structure  causes  the  transverse  branches 
to  split  along  the  main  symmetry  lines  of  the  Brillouin  zone. 
Thus,  modes  from  different  split  TO  and  TA  branches  can 
provide  decay  channels  for  both  the  4 1  (LO)  and  the  E)  (LO) 
modes.  In  the  TA  region,  the  phonon  DOS  presents  a  maxi¬ 
mum  at  100  cm-1,  with  a  significant  contribution  from  the 
lower  TA  branch  along  L-M.  A  secondary  maximum  can  be 
seen  at  —120  cm-1,  with  an  important  contribution  from  the 
nearly  flat  TA  branch  along  L-M.  Similarly,  in  the  TO  region 
the  phonon  DOS  exhibits  a  maximum  at  —455  cm”1  and  a 
secondary  maximum  at  —490  cm”1,  both  of  which  contain 
important  contributions  from  the  split  TO  branches  along 
L-M.  Taking  into  account  the  phonon  DOS,  the  4 1  (LO) 
mode  (574  cm”1)  decays  most  probably  into  pairs  of  TO  and 
TA  modes  with  frequencies  around  455  and  120  cm”1,  re¬ 
spectively,  whereas  the  corresponding  decay  frequencies  for 
the  E,  (LO)  mode  (590  cm”1)  are  clustered  around  490  and 
100  cm-1. 

Assuming  that  the  most  relevant  decay  channel  is  the  gen¬ 
eralized  Ridley  channel  with  the  main  contributions  clustered 
around  oj,  and  w2,  the  temperature-dependent  FWHM  of  the 
phonon  mode  can  be  modeled  by2s 
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FIG.  8.  Raman  spectra  of  the  Aj(LO)  mode  for  temperatures  in 
the  80-560  K  temperature  range,  after  baseline  subtraction.  For  T 
>440  K  the  spectra  were  recorded  with  the  microscope  hot  stage 
and  therefore  the  signal  is  lower  due  to  the  smaller  depth  of  focus  of 
the  setup.  Lower  panel,  fit  to  the  Raman  data  (circles)  used  to  de¬ 
termine  the  baseline  (dotted  line)  to  be  subtracted,  for  the  spectrum 
at  T=  200  K,  where  an  additional  peak  is  already  resolved. 


T(T)  =  r0  +  A[l+n((ol,T)+n{(o2,T)],  (8) 

where  T0  is  a  background  contribution  due  to  impurity 
and/or  defect  scattering  and  isotopic  broadening,  A  is  the 
anharmonic  coefficient,  and  n(a>,T )  is  the  Bose-Einstein  dis¬ 
tribution  function.  Similarly,  the  temperature  dependence  of 
the  frequency  of  the  phonon  mode  can  be  written  as 

cj(T)  =  <o0  +  A0(T)+B[l+n((ohT)  +  n(co2,T)l  (9) 

where  A0(7)  is  the  thermal  expansion  shift  given  by  Eq.  (6) 
and  B  is  the  anharmonic  coefficient. 

As  can  be  seen  in  Figs.  2  and  3,  the  A  ,  (LO)  mode  is  very 
weak  and  emerges  on  a  relatively  strong  second-order  Ra¬ 
man  scattering  background.  Furthermore,  at  low  tempera¬ 
tures  a  peak  at  slightly  higher  frequency  is  resolved.  To  per¬ 
form  a  linewidth  analysis  of  the  A  ,  (LO)  Raman  peak  both 
contributions  must  be  subtracted  from  the  spectra.  In  Fig.  8 
we  show  z(xx)z  Raman  spectra  of  the  A ,  (LO)  mode  for  tem¬ 
peratures  in  the  80-560  K  range  after  background  subtrac¬ 
tion.  To  determine  the  spectra  backgrounds,  a  suitable  com¬ 
bination  of  Lorentzian  line  shapes  plus  a  linear  background 
were  fitted  to  the  spectra,  as  illustrated  in  the  lower  panel  of 
Fig.  8  for  the  spectrum  at  200  K.  The  A  |  (LO)  Raman  peak 
exhibits  a  symmetrical  Lorentzian  line  shape,  which  redshifts 
and  broadens  as  temperature  increases.  For  7’>  400  K  the 
microscope  hot  stage  was  used,  and  hence  the  signal  level  is 
lower  because  of  the  smaller  depth  of  focus  of  the  setup. 
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FIG.  9.  Measured  values  of  frequency  (triangles)  and  linewidth 
(squares)  of  the  A[(LO)  mode  for  temperatures  in  the  80-560  K 
range.  The  dotted  lines  are  fits  of  F(T)  and  oj(T )  given  by  Eqs.  (8) 
and  (9)  to  the  respective  data. 


Measurements  at  T>  560  K  were  discarded  as  the  signal-to- 
noise  ratio  did  not  allow  a  reliable  determination  of  the  peak 
width. 

The  measured  Raman  shifts  and  linewidths  of  the  A ,  (LO) 
mode  are  plotted  in  Fig.  9.  Equations  (8)  and  (9),  with  uq 
=  120  cm-1  and  &>2  =  455  cm-1,  and  T0,  A,  co0,  and  B  as  ad¬ 
justable  parameters,  were  fitted  to  the  respective  experimen¬ 
tal  data.  We  find  r0=1.3  cm-1,  A  =  3.7  cm-1,  (o0 

=  581.4  cm-1,  and  B  =  - 2.1  cm-1. 

A  similar  analysis  was  performed  for  the  £j(LO)  mode. 
Figure  10  displays  the  x(zy)y  Raman  spectra  of  the  /q  (LO) 
mode  for  temperatures  between  80  and  380  K,  after  careful 
subtraction  of  the  background  as  illustrated  in  the  lower 
panel  for  the  320  K  spectrum.  Since  the  Zq(LO)  mode  is 
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FIG.  10.  Raman  spectra  of  the  £j(LO)  mode  for  temperatures  in 
the  80-380  K  temperature  range,  after  baseline  subtraction.  Lower 
panel,  fit  to  the  Raman  data  (circles)  used  to  determine  the  baseline 
(dotted  line)  to  be  subtracted,  for  the  spectrum  at  T=  320  K. 
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FIG.  1 1 .  Measured  values  of  frequency  (triangles)  and  linewidth 
(squares)  of  the  E^LO)  mode  for  temperatures  in  the  80-380  K 
range.  The  dotted  lines  are  fits  of  T(T)  and  co(T)  given  by  Eqs.  (8) 
and  (9)  to  the  respective  data. 

forbidden  in  backscattering  configuration,  the  microscope 
hot  stage  could  not  be  used  and  consequently  this  set  of 
measurements  was  only  possible  on  a  limited  temperature 
range.  Measured  E)  (LO)  linewidths  and  frequencies  as  well 
as  linear-square  fits  of  Eqs.  (8)  and  (9)  to  the  data  with  oj, 
=  100  cm-1  and  co2  =  490  cnT1  are  plotted  in  Fig.  11.  The  best 
fit  to  our  data  is  found  for  r0=2.8  cm-1,  A  =  3.3  cm-1,  w0 
=  595  cm-1,  and  B  =  -l.l  cm-1. 

By  comparing  the  anharmonic  coefficients  A  and  B  for  the 
A  |  (LO)  and  E,  (LO)  modes  of  ZnO  we  can  see  that  the  tem¬ 
perature  dependence  of  these  modes  is  very  similar.  This  is 
not  unexpected,  as  these  modes  lie  close  in  frequency  in  a 
region  where  the  two-phonon  DOS  does  not  change  appre¬ 
ciably.  However,  the  E^LO)  mode  displays  a  larger  line- 
width,  which  implies  a  higher  value  of  T0  and  suggests  that 
the  E^LO)  mode  is  more  strongly  affected  by  impurity 
and/or  defect  scattering  than  the  A  |  (LO)  mode.  As  the  incor¬ 
poration  of  defects  in  ZnO  is  known  to  be  highly  anisotropic, 
this  could  be  related  to  a  possible  anisotropy  in  the  distribu¬ 
tion  of  defects  and  the  fact  that  the  E^LO)  mode  has  an 
in-plane  atomic  motion  whereas  the  atomic  motion  of  the 
A  |  (LO)  mode  takes  place  along  the  c  axis. 

To  sum  up,  in  contrast  with  the  Zs!(lgh  mode  where  anhar¬ 
monic  decay  involves  a  continuum  of  TA+LA  phonons,  the 
temperature  dependence  of  the  A  j(LO)  and  E'](LO)  modes 
can  be  accounted  for  by  a  simple  model  in  which  the  longi¬ 
tudinal  phonon  decays  into  a  pair  of  TO  and  TA  phonons 
whose  frequencies  correspond  to  maxima  of  the  phonon 
DOS. 

3.  Phonon  lifetimes 

The  phonon  lifetime  t  can  be  derived  from  the  Raman 
spectra  via  the  energy-time  uncertainty  relation 

1  A  E 

-  =  — =  27rcT,  (10) 

T  fl 

where  A E  is  the  uncertainty  in  the  energy  of  the  phonon 
mode,  h  is  the  Planck  constant,  and  T  is  the  FWHM  of  the 
Raman  peak  in  units  of  cm-1. 


TABLE  III.  Phonon  lifetimes  of  the  E“gh,  A,(LO),  and  E,(LO) 
modes  of  ZnO  at  80  K  and  room  temperature.  The  characteristic 
decay  time  associated  with  impurities  77  is  estimated  from  the  back¬ 
ground  broadening  parameter  T 0  used  to  fit  the  FWHM  temperature 
dependence. 


Phonon 

mode 

7  = 

80  K 

T= 

300  K 

r0 

(cm  ') 

ti 

(10-12  s) 

1 

(cm-1) 

T 

(10-12  s) 

r 

(cm-1) 

T 

(10-12  s) 

Tjhigli 

3.5 

1.5 

6.0 

0.88 

1.3 

4.1 

AfiLO) 

5.1 

1.0 

9.7 

0.54 

1.3 

4.1 

E,(LO) 

6.6 

0.80 

11.7 

0.45 

2.8 

1.9 

Phonon  lifetime  is  mainly  limited  by  two  mechanisms:  (i) 
anharmonic  decay  of  the  phonon  into  two  or  more  phonons 
so  that  energy  and  momentum  are  conserved,  with  a  charac¬ 
teristic  decay  time  ta  and  (ii)  perturbation  of  the  translational 
symmetry  of  the  crystal  by  the  presence  of  impurities,  de¬ 
fects  and  isotopic  fluctuations,  with  a  characteristic  decay 
time  77.  The  phonon  lifetime  deduced  from  the  Raman  mea¬ 
surements  is  therefore 

-  =  -  +  (11) 

T  tA  Ti 

It  is  difficult  to  separate  the  contribution  of  both  mecha¬ 
nisms,  but  we  can  obtain  an  estimation  of  the  characteristic 
decay  time  associated  with  impurities  from  the  values  of  the 
T0  parameter  provided  by  the  FWHM  temperature- 
dependence  fits.  The  results  are  summarized  in  Table  III.  We 
can  see  that  the  shortest  77  corresponds  to  the  F^LO)  mode. 
The  fact  that  the  F^LO)  mode  has  a  polar  character  allows  a 
Frohlich  interaction  with  charged  impurities  and  defects.32 
As  already  discussed,  owing  to  a  possible  anisotropy  in  the 
distribution  of  defects  such  interaction  could  affect  the 
F^LO)  lifetime  more  severely  than  the  A ,  (LO)  lifetime,  as 
the  E^LO)  mode  has  in-plane  polarization,  whereas  the 
A ,  ( LO)  mode  is  polarized  along  the  c  axis.  This  would  ex¬ 
plain  that  the  E^LO)  linewidth  observed  both  at  80  K  and  at 
RT  is  1.5-2  cm-1  higher  than  the  corresponding  A 1  (LO) 
linewidth. 

The  E3"-"h  lifetime  at  room  temperature  is  in  good  agree¬ 
ment  with  the  value  reported  by  Bergman  et  al.9  No  values 
for  the  A ,  ( LO )  and  E,  (LO)  were  reported  in  that  work.  We 
find  similar  lifetime  values  around  0.5  ps  for  both  modes, 
which  are  close  to  those  reported  for  high-quality  AIN 
crystals.9  The  values  found  for  the  E2lsh,  A ,  (LO),  and 
E^LO)  phonon  lifetimes  corroborate  that  the  anharmonic 
decay  involves  in  all  three  cases  a  three-phonon  process  and 
rules  out  four-phonon  interaction  as  the  main  decay  mecha¬ 
nism,  as  this  would  yield  phonon  lifetimes  2  or  3  orders  of 
magnitude  longer. 

IV.  SUMMARY  AND  CONCLUSIONS 

We  have  carried  out  a  detailed  study  of  the  first-  and 
second-order  Raman  scattering  of  ZnO  for  temperatures 
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ranging  from  80  up  to  750  K.  Second-order  Raman  spectra 
obtained  in  different  scattering  geometries  have  been  dis¬ 
cussed  in  terms  of  DFT  phonon  dispersion  relations  and 
symmetry  analysis.  Characteristic  second-order  Raman 
peaks  have  been  assigned  to  overtones  and  combinations  of 
acoustic  and  optical  phonon  branches.  The  temperature  de¬ 
pendence  of  the  333  cm-1  Raman  peak,  which  is  often  attrib¬ 
uted  to  acoustic  overtones  at  M,  has  been  shown  to  corre¬ 
spond  to  a  difference  mode,  and  the  mode  has  been  assigned 
to  £,2Igh-£'1°w.  Several  sharp  peaks  show  up  at  low  tempera¬ 
ture,  which  are  tentatively  related  to  intrinsic  lattice  defects. 
Most  notably,  a  peak  is  resolved  at  589  cm-1  in  the  z(xx)z 
spectra  at  80  K,  very  close  to  the  Ztj(LO)  frequency.  The 
possibility  of  a  leaky  Ztj(LO)  mode  can  be  ruled  out  since  the 
frequency  of  the  Zs^LO)  mode  at  80  K  was  determined  to  be 
593  cm-1  in  the  allowed  x(zy)y  configuration. 

The  anharmonic  broadening  of  the  E 2lgh,  A  |  (LO),  and 
Zij(LO)  phonon  modes  was  studied  as  a  function  of  tempera¬ 
ture.  The  £*lgh  mode  of  ZnO  exhibits  a  visibly  asymmetric 
line  shape,  which  could  be  successfully  explained  in  terms  of 
resonant  anharmonic  interaction  of  the  Z?,lgh  mode  with  a 
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band  of  combined  transverse  and  longitudinal  acoustic 
modes,  as  the  steep  variation  of  the  two-phonon  DOS  around 
the  Zs,lgh  frequency  leads  to  a  distorted  phonon  line  shape.  In 
contrast,  the  two-phonon  DOS  is  rather  flat  in  the  A  |  (  LO) 
and  Zs^LO)  frequency  region  and  the  temperature  depen¬ 
dence  of  these  modes  can  be  adequately  described  with  a 
simple  anharmonic  decay  model  via  a  dominating  Ridley 
channel  involving  TO  and  TA  modes  whose  frequencies  are 
clustered  around  maxima  of  the  DOS.  The  lifetimes  of  these 
modes  are  determined  from  the  Raman  spectra  and  are  found 
to  be  around  0.5  ps  for  both  longitudinal  modes  at  room 
temperature. 
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